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Abstract. We study the long term evolution of the relationship between the root mean square (rms) variability and flux (the 
"rms-flux relation") for the black hole Cygnus X-1 as monitore d from 1996 to 2003 with the Rossi X-ray Timing Explorer 
(RXTE). We confirm earlier results by Uttlev & McHardv ( 2001 ) of a linear relationship between rms and flux in the hard state 
on time scales > 5 s reflecting in its slope the fractional rms variability. We demonstrate the perpetuation of the linear rms-flux 
relation in the soft and the intermediate state. The existence of a non-zero intercept in the linear rms-flux relation argues for two 
lightcurve components, for example, one variable and one non-variable component, or a possible constant rms component. The 
relationship between these two hypothesized components can be described by a fundamental dependence of slope and intercept 
at time scales < lOksec with long term averages of the flux. 
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1. Introduction 

Galactic Black holes (BHs) are generally found in two ma- 
jor states which are mainly characterized by t heir X-ray and 
radio spectral and timing properties (see, e.g.. iNowakl I2003. 
and references therein). The "hard state", occuring at lumi- 
nosities L < 0.05 LEdd, where Lgdd is the Eddington luminos- 
ity, is characterized by ~30% root mean square (rms) vari- 
ability, a hard X-ray spectrum that can be roughly described 
by an exponentially cutoff power law, and radio emission. 
The hard state energy spectrum can physically be modelled 
by thermal Comptonization of cool seed photons in a hot 
electron plasma (Thorne & Price, 1975; Shapiro et al., 1976; 
ISunvaev & Trumped. 1 1979llDove et allll998l) . At higher lumi- 
nosities, BHs are seen in the "soft state" with a thermal X-ray 
spectrum with a characteristic temperature of at most a few 
keV, X-ray variability of only a few percent rms, and no ob- 
served radio emission frendedkood) . The thermal X-ray spec- 
trum in the soft state is generally associat ed with a geometri- 
cally t hin, optically thick accretion di sk (Novikov & Thorne, 
1 19731 IShakura & Sunvaevt Il973t iGierlinski et all Il998h . In 
addition to these two canonical states, further states have been 
identified which are characterized either by an even larger lu- 
minosity than in the soft state (the "very high state"), or by 
variability and X-ray spectral properties which are mostly in- 
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termediate between the hard and the soft state (the "intermedi- 
ate state": iBelloni et all fl996). 

The analysis of the short term X-ray variability of BHs is 
especially well suited for gaining better insight into the physi- 
cal processes at work near the compact object. In order to an- 
alyze this behavior in a systematic way, in 1997 we initiated 
a long term monitoring campaign of the canonical hard state 
BH Cygnus X-1. During the campaign, which is still ongoing, 
Cyg X-1 is observed at ~7d to ~14d intervals simultaneously 
with the Rossi X-ray Timing Explorer (RXTE), the Ryle radio 
telescope, and several optical instruments. These observations 
enable us to study the evolution of the observational properties 
of Cyg X-1 (e.g., X-ray spectrum and timing parameters, radio 
flux) and thus the source history on time scales longer than the 
orbital time scale of ~ 5.6 days. 

Earlier results from our ca mpaign have been reported by 
IPottschmidt et alJ ilOOdh and bv lPottschmidt et all J20031 here- 
after paper I), where we studied the evolution of the Fourier- 
frequency dependent X-ray time lags an d of the power spectral 
density (PSD), respectively. Following Nowak (2000), in pa- 
per I we were able to show that for a very wide range of source 
parameters, the PSD of Cyg X-1 is well described by the su- 
perposition of broad Lorentzian functions. Such a description 
is an alternative to the more c ommon desc ription of the PSD in 
terms of broken power laws jvanderKlisU 19951) or shot-noise 
profiles llLochner et all ll 99 lL and references therein). During 
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changes of the source into the "intermediate state" or the "soft 
state", the characteristic frequency of these components shifts 
towards higher frequencies in a systematic manner, some of 
the Lorentzians change their strength or vanish, and the X-ray 
time lag increases in the frequency b and dominated by those 
Lorentzians which become stronger (Pott schmidt et al j. [2000, 
paper I). Once the "soft state" is reached, the X-ray time lags 
are again similar to those seen during the hard state. 

In paper I we interpreted these results in terms of 
models invoking damped oscillat ions in the accretion disk 
and/ or the Comptonizing plasma ( | di Matteo & Psaltis. 1999; 
Psaltis & Norman, 200 lj IChurazov et all [200 1 ; No wak et all 
2002). Because of the complicated physics of the system "ac- 
cretion disk - accretion disk corona - radio outflow", none 
of these models are fully self-consistent and thus none of the 
models describes all observed properties of Cyg X-l. It is our 
hope that additional empirical data for the hard state behavior 
might help to further constrain these different emission models 
ofBHs. 

Recently, lUttlev & McHardvl d200lb found a remarkable 
linear relationship between the rms amplitude of broadband 
noise variability and flux in the X-ray lightcurves of neu- 
tron star and BH X-ray binaries and active galactic nuclei 
(AGN). This linear relation between flux and rms variability 
has an offset on the flux axis, leading to the suggestion that 
the lightcurves of X-ray binaries are made from at least two 
components: one component with a linear dependence of rms 
variability on flux, and one component which contributes a con- 
stant rms to the lightcurve or does not vary at all. The fact that 
the same rms-flux relation was found in the lightcurves of ac- 
tive galactic nuclei suggests that this behavior is intrinsic to 
compact accreting systems. 

In the case of X-ray binaries, the rms-flux relation has 
until now only been investigated for single observations of 
SAX J1808.4-3658 and for a few hard state observations of 
Cyg X-l. It is important to discover whether this relation ap- 
plies generally for many observations of a given object and 
across a range of states. Therefore in this paper we investigate 
the evolution and behavior of the rms-flux relation in Cyg X- 
1 using the RXTE monitoring data. Preliminary results have 
already been presented elsewhere Ideissner et allE003t) . 

The remainder of this paper is structured as follows. In 
Sect. 13 we describe the data and the computation of the rms 
variability. We then present our results on the rms-flux relation 
(Sect. [3}. We show that the rms-flux relation is valid through- 
out all spectral states, for all energies and on all time scales. 
It is demonstrated how mean flux, (F), and the relation fit pa- 
rameters slope, k, and intercept on the flux axis, C, map out 
a fundamental plane in the hard state. After analyzing spectral 
dependences, we examine the rms-flux relation on short time 
scales by calculating the variability of the rms variability. We 
end the paper with a discussion of the physical interpretations, 
and explore the meaning of the long term rms-flux relation as 
well as of the hard state k-C-{F) fundamental plane (Sect.|4}. 



Table 1. Definition of the energy bands and data modes. 
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2. Observations and data analysis 

2.1. Data extraction 

Most of the data analyzed in this paper are the lightcurves that 
we already used in paper I of this series. We therefore only give 
a brief summary of the RXTE campaign and the data extraction 
issues and refer to paper I for the details. 

The data analyzed here spans the time from 1996 until 
early 2003, covering our RXTE monitoring programs since 
RXTE Announcement of Opportunity 3 (A03, 1998) through 
RXTE's A06 (2002/2003). Additionally, other public RXTE 
data from RXTE AOl and A07 were used. Typical exposure 
times were 3 ksec in 1998, and lOksec since then. After screen- 
ing the data for episodes of increased background (see paper I), 
we extracted lightcurves with a resolution of 2~ 6 s (-16 ms; 
A03) or 2 -8 s (~4ms; other data) using data fro m the RXTE 
Proportional Counter Array ( Jahod a et all 1996) and the stan- 
dard RXTE data analysis software, HEAS0FT, Version 5.0. We 
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extracted lightcurves in 5 different energy bands (Table^. Due 
to integer overflows in the binned data modes, high soft X- 
ray count rates can lead to a distortion of the rms-fiux relation 
for the lowest energy band (see Appendix A for further dis- 
cussion). We therefore use only data from energy bands 2-5 
which are not affected. As the total number of operational pro- 
portional counter units (PCUs) was often not constant during 
an observation, separate lightcurves were generated for each of 
the different PCU combinations. These lightcurves were ana- 
lyzed separately. To facilitate comparisons, unless noted other- 
wise we normalize all data to one PCU. 

Tabled shows that the energy limits of the high-resolution 
lightcurves, which were used for the computation of the rms- 
flux relation, in parts differ significantly. In order to enable 
comparisons, we normalize all observations to the energy lim- 
its of epoch 4 (P40099-06/28, P50110, P60090, P70414). We 
calculate a flux correction factor to the epoch 3 observations by 
integrating the count rates of the corresponding spectrum chan- 
nels of each observation. This factor has to be applied also to 
all values that are dependent on flux, e.g. C, cr, etc. 

2.2. Computation of the rms variability 

To determine the rms variability as a function of source flux, 
we first cut the observed lightcurve into segments of 1 s dura- 
tion for which the source count rate is determined. These seg- 
ments are then assigned to their respective flux bins. In general, 
40 linearly spaced flux bins were chosen. In the final compu- 
tation, only flux bins containing at least 20 segments were re- 
tained. We note that short time segments of 1 s duration are 
necessary to study the characteristic variability of the source 
which is dominated by frequencies above 10 _I Hz. Thus us- 
ing a 1 s segment size allows us to measure rms-fiux relations 
which sample a broad range of fluxes, while also measuring 
rms with high signal-to-noise, since the minimum frequency 
of variability sampled within each segment is well below the 
frequency where photon counting noise dominates. 

Next for each of the lightcurve segments in each of the flux 
bins, the rms variability is determined for the Fourier frequency 
interval of interest. We first compu te the PS Ds of the ind ividual 
lightcurve segments (iNowak et allll99 9: Nowak, 2000), using 
the PSD normalization where the integrated PSD equals the 
squared fractional rms variability (i.e., the fractional variance) 
of th e lightcurve ( Bell oni & Hasingen. Il990t iMivamoto et all 
Il992l) . We then average all the PSDs obtained for segments in 
the same flux bin, and bin the averaged PSD over the frequency 
range of interest (for this paper we use the range from 1 Hz 
to 32 Hz) to yield the average power density in that frequency 
range, (P). Then, the absolute rms variability of the source in 
the frequency band of interest, cr, is obtained according to 



cr= [«P> - c Poisson ) • Af] l/2 F 



(1) 



where cp i SS on is the Poisson noise level due to photon counting 
statistics, Af is the width of the frequency band of interest, and 
where F is the source count rate of the flux bin (normalized to 
one PCU). Note that in the comparatively low frequency range 
considered here, the detector deadtime does not need to be 
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Fig. 1. Relation between the PCA count rate, F, and the rms 
variability, cr, for the observation of 1998 July 02, showing the 
typical hard state characteristics. The data are extracted from 
energy band 2 (4.1-6.2keV). The dotted line is the best-fitting 
linear model to the data (see text), the error bars are at the lcr 
level. 



taken into account when determining the Poisson noise level. 
The statistical uncertainty of the average PSD value, A(P), 
is de termined from the periodogram statistics ll van der Klis. 
1989) 

<"> (2, 



A<P> 



Vmw 



where M is the number of segments used in determining the av- 
erage rms-value and W is the number of Fourier frequencies in 
the frequency range of interest. The uncertainty of cr, A(cr), is 
then computed using the standard error propagation formulae. 
Note that, since we measure cr from the average PS D in each 
flux bin, our method differs from that of Uttlev & McHardv 
(2001), who first calculated a Poisson noise-subtracted rms for 
each individual segment before binn ing up to obtain cr in each 
flux bin. Uttlev & McHardv (2001) computed the uncertainty 
of cr by calculating standard uncertainties from the scatter of 
the rms values in each flux bin. This method has the disadvan- 
tage that, when the variance due to the source is small com- 
pared to the variance due to Poisson noise, subtracting the ex- 
pected Poisson noise level can sometimes lead to negative vari- 
ances due to the intrinsic variations in the realized noise level. 
Our revised method avoids this problem by averaging the seg- 
ment PSDs before noise is subtracted, thus reducing the intrin- 
sic variations in the noise level. 

2.3. rms versus flux 

Fig-fflshows the rms versus flux relationship for a typical hard 
state observation, taken on 1998 July 02, in energy band 2 (4. 1- 
6.2 keV). To describe th is linear relationship, we model it with 
a function of the form JUttlev & McHardvthoOll) 



cr = k(F-C) 



(3) 
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where the slope, k, and intercept, C, are the fit parameters. 
This linear model fits the general shape of the rms-flux re- 
lation of the example reasonably well, although the large x 2 
value of 102.1 for 22 degrees of freedom implies that there 
is some weak intrinsic scatter. The best-fitting model param- 
eters for the example of Fig. ^ are k = 0.30 ± 0.01 and 
C = 67 ± 5 counts s _I PCLT 1 (unless otherwise noted, uncer- 
tainties are at the 90% confidence level for two interesting pa- 
rameters). 

In order to interpret the parameters of Eq. it is useful to 
use a toy model in which the observed source flux, F obs(t), is 
written as the sum of two components, 

F ohs (t) = Fconst + ^var(0 (4) 

where F const is a component of the lightcurve which is as- 
sumed to be non-varying on the time scales considered here, 
and where F val (t) is a component which is variable and thus 
responsible for the observed rms variability. In the toy model 
of Eq. © the interpretation of the slope k of the rms-flux re- 
lationship is straightforward: it is the fractional rms of F var (f). 
Furthermore, the intercept of the rms-flux relation on the flux 
axis, C, is the count rate of F const . It is obtained by extrapolating 
F var — > in Eq. ®. 

While the decomposition of Eq. ® is the most straight- 
forward, it is not fully unique. As we will show later, for the 
majority of our observations C > 0, however, there are several 
observations in which C < 0. In these cases the above interpre- 
tation clearly breaks down (there are no negative fluxes) and it 
makes more sense to write Eq. (0 as 

cr = kF + o-o (5) 

In this formulation of the rms-flux relationship we have intro- 
duced ctq, an excess rms present in the lightcurve, which is 
positive whenever C < 0. Note that for ctq > Eq. (0 pre- 
dicts cr — > ctq for F — > 0, i.e., the extrapolation of the rms-flux 
relationship to flux zero predicts that there is still variability 
present, which is clearly unphysical. We do not consider this 
behavior a problem, however, as the rms-flux relationship is 
measured only for F » 0. We deem it reasonable to assume 
that for F - other variability processes which do not obey the 
linear rms-flux relationship will become important and force 
cr — > for F — > 0. Finally, we note that for ctq < the ob- 
served rms is reduced with respect to a linear rms-flux rela- 
tionship with (To = 0. This is possible, e.g., if a constant flux 
component, F const , is present in the lightcurve (the majority of 
our fits imply reductions of about 25% in rms with respect to 
cro = 0). 

Observationally, degeneracies in the fit prevent us from dis- 
tinguishing between the two interpretations of the rms-flux re- 
lationship. Since most of the observations of Cyg X-l pre- 
sented in the following have C > 0, and to enable the com- 
parison with other analyses, we have chosen to use Eq. Q as 
our fit function. We ask the reader, however, to keep the second 
interpretation of the rms-flux relationship in mind. 

The best test of linearity is a \ 2 test of whether a straight 
line fits the data. On the other hand, it is usually very diffi- 
cult to fit a simple model to variability properties because of 



the scatter in the data points. It may be misleading to discard 
an observation as being 'non-linear' just because a straight line 
does not supply a good^ 2 fit. We found that Kendall's rank cor- 
relation coefficient, which tests whether a correlation is mono- 
tonic, serves as a good indicator of linearity in the sense that the 
bulk of points follows a straight trend, while allowing for some 
points, predominantly at high fluxes (cf. Fig.^, to diverge. 

Therefore we choose to characterize the goodness of the 
linear correlation between N data points of cr and F using 
Kendall's rank correlation coefficient, t, defined by 



N(N - l)/2 

Here S is the sum of scores which are determined if we assign 
to each of the N(N - 1) /2 possible pairs [(cr-JFi), (Oj/Fj)] of 
the sample (cn/Fi), . . . , (ctn/Fn) a score of +1 or -1 depend- 
ing on whether their ranks are in the same or der or in the o ppo- 
site order on the cr and F axes, respectively (Keeping, 1962). 

The correlation in Fig.[0 e.g., has r = 0.95 (the range of the 
correlation coefficient is -1 < t < 1, with a perfect correlation 
being indicated by t = 1). For the remainder of this paper we 
set a minimum threshold of T m ; n = 0.9 below which we do not 
accept the hypothesis of a linear rms-flux relationship. Given 
this criterion, the fit with this simple linear function works for 
the majority of all observations (201 out of 220). With a few 
exceptions, the rejected observations with t < r m m generally 
have short exposure times and therefore bad counting statis- 
tics, resulting in a scattered rms-flux relation. In the following 
discussions of the behavior of k and C we will not include these 
observations (but see Sect. 13. 21 . 

Even though the^ 2 fits for linearity are not always so good, 
the maximum deviations from the best fit straight line are gen- 
erally less than 15%. We therefore think that it is acceptable to 
test linearity by r and use the fit parameters k and C for correla- 
tions with other variables, although these correlations may con- 
tain some approximate only data points due to fitting a model 
which does not exactly describe the data. 

3. Results 

3.1. Evolution of slope and intercept 

Figs. |2^-c show Kendall's t, the mean RXTE PCA count 
rate, (F), for energ y band 2 (-4-6 keV ), and the RXTE All 
Sky Monitor (ASM: lLevine et alll996l) count rate for the time 
from 1996 to 2003. In Figs. EJj and c the soft state that took 
place 1996 March-September and the "soft state" of 2001 
October-2002 October are reflected by increased count rates. 
The fact that the failed transitions and the soft state behavior of 
2001/2002 look similar in C together with the observation that 
the radio emission is not always quenched in the 2001/2002 
soft state (G.G. Pooley, private communication), allows one 
to speculate that the 2001/2002 s oft state s hows intermediate 
state phases (see, e.g., Zdzia5kHII][2002) 

. Given the uncer- 
tain nature of the 2001/2002 soft state we must be careful not 
to assume that it is a canonical soft state or equal to the 1996 
soft state without further proof (Wilms et al. 2003, in prep.). 
Between the 1996 soft state and the 2001/2002 interval, Cyg X- 
1 was usually found in the hard state, infrequently interrupted 
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Fig. 2. a) Kendall's r, the long-dashed line indicates r m j n , observations with t < T m [ n were not included in the further analysis, 
b) Mean RXTE PCA count rate for en ergy band 2 (~4-6keV) . c) Mean RXTE ASM count rate, rebinned to a resolution of 
5.6days, the orbital period of Cyg X-l ( Br ocksopp et all 11999). The four most prominent intermediate states are indicated by 
vertical dotted lines. 
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Fig.3. Fit parameters and Cyg X-l system properties for energy band 2 (~4-6keV). a) Slope, and b) x-axis intercept of the 
rms-fiux-relationship, c) mean rms value, (cr), d) mean RXTE PCA count rate for energy band 2. Error bars given for slope 
and intercept are uncertainties at the 90% confidence level for two interesting parameters. Triangles indicate observations which 
cannot be described by a linear rms-fiux relation (i.e., r < T m ; n ). 



by the "failed state transitions" described in the introduction 
and paper I (see also Sect. 13. 21 . 

Figs.[5k-d show the fit parameters k, C, mean absolute rms 
variability, (cr), and mean RXTE PCA count rate, (F). The 
most apparent correlation is between (cr) and (F). To examine 



this correlation, we display the rms-fiux relation of the long 
term variations in cr in Fig. [4] A tight linear correlation between 
(cr) and (F) can be seen in the hard state and, with a slightly 
lower slope, in the soft state data. It seems that the linear rms- 
fiux relation observed on short time scales also applies on much 
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Fig. 4. Mean rms variability, (cr), as a function of mean PCA 
count rate, (F), The hard state observations line up to follow 
a linear rms-flux relationship. There is also a linear rms-flux 
relationship for the soft state observations, but with a slightly 
smaller slope. 



longer time scales. A least x 1 fit ix 1 - 15589 for 161 dof) of a 
linear model to the hard state provides k = 0.2230 + 0.0003 and 
an intercept of C — 6.8 + 0.2 counts s _1 PCU -1 . These values 
of k and C are lower than the values typically observed in the 
short term rms-flux relation. The corresponding fit for the soft 
state data (x 1 = 30660 for 50 dof) provides a slope of k = 
0.1867 ± 0.0004 and an intercept at the flux axis of C = 16.8 ± 

0. 6 counts s PCU -1 . The values for the failed state transitions 
are associated with either the hard state or the soft state data. 

We next consider the long term behavior of k (Fig. [3Jt). 
During the soft states of 1996 and 2001/2002, k takes a sig- 
nificantly lower value than during the hard state. During the 
last three "failed state transitions" before the 2001/2002 soft 
state k decreases to the soft state level. The values of the inter- 
cept show different levels, too: while in the 1996 soft state C is 
larger than or equal to the hard state, C in the 2001/2002 "soft 
state" is smaller than in the hard state. This may be connected 
to the special soft state character of the 2001/2002 interval. The 
detailed characteristics of the rms-flux relation in the different 
states will be discussed in the following sections. 

One of the main results of paper I was that in 1998 May a 
change of the general long term behavior of Cyg X-l from a 
"quiet hard state" to a "flaring hard state" took place that coin- 
cided with a change in the PSD shape, resulting in a decrease 
of the relative rms amplitude from an average of 36 ± 1 % to 
29 ± 1%, for the total 2-13.1 keV power spectrum (see paper I, 
Fig. 3d). This change is seen in Fig.[3]as a rise of k and C, while 
the mean rms value, (cr), stays on a constant level. This appar- 
ent discrepancy can be explained by taking into account that the 
total rms is measured from ~ 2 x 10~ 3 Hz to 32 Hz in paper I, 

1. e., over a much broader frequency range than here. As we 
show below (Sect. 14. 2l . the drop in the global rms in 1998 May 
is due to a change of the PSD outside of the frequency range 
covered by our measurements of the rms-flux relationship. 
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Fig. 5. Contour plot at the 90% confidence level showing the 
two fit parameters slope and intercept at the flux axis. 

To illustrate the correlation between k and C, Fig.|5]shows 
the 90% confidence contours for these parameters. The points 
in the £-C-plane fall into two separate regions corresponding to 
hard and soft states. A possible correlation between k and C can 
be seen within the soft state, but a much stronger correlation 
between k and C is present within the hard state. Interestingly, 
the k-C values in the hard state fall into a distinctive, fan-like 
shape, which suggests that a more fundamental relation may 
underly the k-C correlation in that state. The "fan" appears to 
converge at the values of k and C measured for the long term 
rms-flux relation of the hard state (Fig.©, suggesting that the 
pattern of k and C observed in the hard state may be related to 
the long term rms behavior. 

In fact it is easy to see how a linear long term rms-flux re- 
lation can produce the observed fan-like shape of the k-C cor- 
relation in the low state. First, consider that the mean cr of the 
z'th observation, (cr,), is related to the mean flux, (Ft), by 



<try> = ki((Fi) - Q 



(7) 



where fc, and C, are the k and C values determined from the 
rms-flux relation of that observation. According to Fig.|4] (cr,) 
is also given by 



<CT,) = k long ((Fi) - Cong) 



(8) 



where k\ ong and Cong are the k and C values of the long term 
rms-flux relation (k\ ong = 0.2230, C ong = 6.8). From Eqs. 
and (jSJl it is easy to see that 



ki = 



and 



c 



^long {(J*i) ClongJ 

<F,-> - C 



k t (F i) - k iong (<F,-> - Cong) 



(9) 



(10) 



The distribution of the C values is thus caused by the 
requirement to maintain the rms-flux correlation on all time 
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Fig. 6. For hard state observations, k, versus ((Fj) - 
Ciong)/ {{Fi) - Ci) gives a linear relation y — bx with slope 
b = 0.2312. The error bars on k are the formal lcr error, the 
horizontal error bars are calculated from propagating formal 
lcr errors of Fj and C,. 

scales. The "fan like" distribution of the individual values is 
the result of the requirement that for k t — > k\ ong , Q — > Ci ong . 
Confirming earlier results, we find that k-, is not correlated with 
flux. In order to maintain the long term rms-flux correlation, 
this results in a corresponding scatter of the values of C, , form- 
ing the "fan". We can test this directly by plotting ki versus 
((Ft) - Ciong)/ ((Fi) - Ci) (Fig.|6j, which gives a linear relation 
Or 2 = 308 for 153 dof) with slope 0.2312 ± 0.0004 £ long ), 
as expected from Eq. We conclude that the various values 
of k, C and mean flux (F) map out a "fundamental plane" in 
the low state, governed by the requirement that the long term 
linear rms-flux relation is maintained. We will examine the con- 
sequences of this fact in detail in Sect. @] 

Now that we have described the general evolution of the 
parameters of the rms-flux relationship, we turn to describing 
the individual relations in greater detail. 

3.2. The rms-flux relation in the hard and soft state 
and during "failed state transitions" 

The high values of Kendall's r in Fig. |2^ prove that a linear 
rms-flux relation is valid throughout all states, the hard state 
and the soft state and the "failed state transitions", with rela- 
tively few outliers. Note, however, that a number of soft state 
observations, predominantly during the atypical 2001/2002 
"soft state", have rms-flux relations with a negative flux off- 
set C (see Fig. |3J>). This negative offset cannot be explained 
if, for these cases, C simply represents a constant component 
to the lightcurve, but can easily be explained if C represents a 
component with constant rms. In Fig.^-c we give examples 
for the overall rms-flux linearity. 

There are a small number of observations which cannot be 
described with a linear model. These deviations from a lin- 
ear relation have two principal reasons: (1) As discussed in 
Sect. 12. 31 many observations with r < r m i n simply have a short 



observation time, so that the number of segments per bin is rela- 
tively low which produces the strong scatter, resulting in a poor 
correlation. (2) In contrast, the "wavy" rms-flux relation shown 
in Fig.lSJ) results from short term variations in the parameters C 
and/or k of a linear rms-flux relation. If a single observation (on 
time scales of a few 100 s) consists of several parts, each with 
different average flux (see Fig. [SJi), then the superposition of 
the linear rms-flux relations of all parts will result in a "wavy" 
rms-flux relation, as each part displays a slightly different C 
and/or k value. 

To describe the flaring events when Cyg X-l attempts the 
change from the common hard state into the soft state, without 
reaching the soft state, we coine d the term "failed state tran- 
sitions" ( Pottschmi dt et alJ l2000: paper I). Depending on how 
far its state evolves, during these failed transitio ns Cyg X- 1 ca n 
reach the intermediate state as defined bv lBelloni et alJ fi*996). 
or its behavior can remain hard state-like during these flares. 

We examined four "failed state transitions" that were cov- 
ered by RXTE observations in detail: 1998 July 15, 1999 
December 05 (extracted as two separate lightcurve parts), 2000 
November 03, and 2001 J anuary 29 (probably a short soft 
state, see ICui et all l2002h . The values of k and C of the 
first failed state transition on 1998 July 15 are not different 
from the surrounding hard state observations. During the lat- 
ter three failed state transitions, however, the values of k and 
C change to the soft state level. For one exemplary observation 
on 2000 November 03, this behavior is displayed in Fig. 05. 
Comparing the failed state transition with the hard state obser- 
vation immediately before (Fig. Et), the rms-flux relation of 
2000 November 03 changes to a lower slope and a small, neg- 
ative intercept C during the "failed state transition" itself. The 
PSD of this observation shows that this change in the rms-flux 
relation is accompanied by a decreasing Lorentzian component 
L\ and an increasing Lorentzian component L3 (Fig.01-e). 

3.3. Spectral dependence of the rms-flux relation 

We confirm a correlation between the photon spectral index, T, 
i.e., the hardness of the spectrum, and the mean rms, (o~), of 
the examined observations. In the hard state, the photon spec- 
trum of each observation can roughly be described as the sum 
of a power law spectrum E~ T with photon index T and a multi- 
temperature disk black body ( Mitsuda et al ., 1984). The corre- 
lation of Fig.[5k takes into account only hard state and "failed 
state transition" observations and shows that (cr) increases as 
the spectrum grows softer. As rms variability and flux are lin- 
early correlated in the hard state, this correlation between F 
and (cr) is equivalent to a correlation between T and flux. A 
similar result with RXTE ASM data of Cyg X-l has been 
recently shown by Zdziarski et al 1 d2002h . who found a very 
strong hardness-flux anticorrelation in the hard state. 

We note that we are using the definition of paper I, using the 
term "failed state transitions" for those observations which ex- 
hibit a significantly increased time lag compared to its typical 
value, since the X-ray time lag seems to be a more sensitive in- 
dicator for state changes than the X-ray spectrum. As has been 
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already been noted in Sect. l3.2l there are only four correspond- 
ing observations. 

The filled diamond at T ~ 2 in Fig.|9]belongs to the "failed 
state transition" of 1998 July 15. The fact that this data point 
is situated within the value range of hard state observations in 



Figs. E] and ^confirms that this "failed state transition" has to 
be treated separately and classified as being close to the hard 
state, as already mentioned in Sect. 13.21 Allowing for the split 
in values of k and C, respectively, between the hard state and 
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Fig. 9. a) Mean rms variability, {&), b) slope, k, and c) flux in- 
tercept, C, as a function of photon spectral index, F, for failed 
state transitions and the hard state observations. Diamonds des- 
ignate "failed state transitions", crosses are for hard state obser- 
vations. Error bars for k and C are given at the 90% confidence 
level for two interesting parameters. 



failed state transition data, any correlation between the photon 
index and k, or C, cannot easily be detected. 

So far we concentrated on the rms-flux relation in the sec- 
ond lowest energy band, i.e., energy band 2 in the range ~4— 
6keV. This approach is justified as the linear rms-flux relation 
is found to be valid for all energy bands considered here (see, 
e.g., Fig. I10i. such that the general behavior described for en- 
ergy band 2 is also seen in the other bands. Even the rms-flux 
relation in energy band 3 which contains the iron K line does 
not show a different behavio r, in agreement with the results of 
iMaccarone & CoppH J2002al) that the iron line tracks the con- 
tinuum at least in the soft and transition states. Nevertheless, 
there are subtle differences in the shape of the rms-flux corre- 
lation. 

Comparing the values of k and C for the different energy 
bands, our energy resolved results for the rms-flux correlation 
show that k, which represents the fractional rms variability of 
the time variable component of the emission in the 1-32 Hz 
band, generally decreases with energy when the source is in 
the hard state (Fig.llH . in agreement with earlier results (e.g., 
iNowak et allll998l) . The same is true for C (see Fig. ll2l >. For 
the highest energy band 5, our analysis reveals that k and C 
only slightly correlate with energy band 2 (Figs.fTTV and!12b). 
This result confirms earlier work claiming only a 'loose cou- 
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Fig. 10. Shape of a)-d) the rms-flux relation and e)-h) the cor- 
responding PSD for the observation of 1999 May 21 for energy 
bands 2-5, respectively. 



pling' b etween the source variability at low and at high energies 
(Zv ckil 120031 biurazov et all 12001 1 IMaccarone et all [2000; 
Gilf anov et allll999l) . In terms of Compton corona models, the 
fraction of photons from the accretion disk significantly de- 
creases when proceeding from energy band 2 to energy band 5. 
Therefore we take our result to confirm earlier claims that the 
variability properties of Cyg X-l are driven by coronal fluctu- 
ations and not by chang es in the soft photon input to the puta- 
tive C ompton corona (Macc arone et al ..2000: IChurazovetall 
l200ll and references therein). 

To examine the overall spectral shapes of the two variability 
components characterized by C and k, respectively, we plotted 
these values for the energy bands 2-5, normalising C by (F) 
to account for the instrument response. Fig.ll3lshows two ex- 
amples with different photon spectral index, F. For most obser- 
vations with T lower than ~ 2.0, which means for the majority 
of hard state observations, the spectral shapes of C and of k are 
similar, but the spectral shape of C tends to be somewhat flatter: 
both C and k show a soft energy spectrum, i.e., higher values at 
low energies and lower values at higher energies (Figs.[T3t-b). 
When we proceed to higher F, observations begin to display a 
flat or even tilted spectrum, i.e. lower values at low energies 
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and higher values at higher energies, more pronounced for C 
than for k (Figs. 113b— d'). The change in the shape of the nor- 
malised C (cf. Figs. I13h and c) is caused by the combination 
of an intrinsic change in the spectral shape of the unnormalised 
C and a change in the spectral shape of (F). Suggestions that 
the two variability components C and k, representing a constant 
and a variable rms component, will have basically different en- 
ergy shapes cannot be confirmed, but to a certain degree in the 
low state the spectral shape of C is flatter and harder than the 
spectral shape of k. 



3.4. The rms-flux relation on short time scales 

The fact that <x tracks the flux implies that the "lightcurve" of 
the rms values is similar to the flux lightcurve. This fact can be 
used to probe how the rms responds to flux changes on short 
time scales (see also Uttley et al. 2003, in prep.). For example, 
if the rms does indeed track variations on all time scales, the 
PSD of the "rms lightcurve" should look similar to the PSD of 
the conventional "flux lightcurve". In contrast, if the rms only 
tracks variations on time scales longer than, say, 10 s, the PSD 
of the rms lightcurve should be sharply cut off at 0. 1 Hz. In or- 
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Fig. 13. Spectral shapes of the two variability components char- 
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(F). a)-b) show a typical hard state observation (2001 April 
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Fig. 14. PSDs of a) F lightcurve and b) cr lightcurve for 
the hard state observation on 1999 July 05. The original flux 
lightcurve with time resolution of At = 2~ 8 s has been di- 
vided into 0.25 s segments, and the frequency range v = 4- 
128 Hz was chosen for calculation of the rms variability of 
each segment. To compute PSDs from the flux and the rms 
lightcurves - consequently both with the time resolution of the 
segment length At = 0.25 s - segments of 1024 time elements 
(= 1024 x 0.25 s = 256 s) were used, thus resulting in a fre- 
quency range from v m i n = 0.004 Hz to v max = 2 Hz. We re- 
binned this range to 84 logarithmically spaced bins, c) Ratio 
of the F lightcurve PSD to the cr lightcurve PSD in order to 
quantify the similarity of the two PSDs. 



der to test this issue we created 0.25 s lightcurve segments and 
calculated F and cr for each of these segments to make flux 
and rms lightcurves. From these lightcurves, we calculated the 
PSDs (shown for an arbitrarily chosen observation in Fig. I14i. 
In this case, we do not correct the rms lightcurve for Poisson 
noise, since the stochastic nature of the noise process means 
that such a correction is only possible when averaging over 
comparably long time intervals. Note, however, that the effect 
of the Poisson noise and intrinsic noise contributions to the rms 
lightcurve is to add power to the rms PSD at high frequencies, 
not to cause spurious high-frequency cut-offs or breaks. 

Over the frequency range ~0. 01-0. 2 Hz the rebinned PSDs 
show the same features as indicated by the flat ratio in Fig. 114b. 
Therefore it is likely that the rms-flux relation is fulfilled on 
time scales from 5 s to 100 s. Above 0.2 Hz, however, Poisson 
noise starts to dominate the rms PSD such that the shape of the 
rms and flux PSDs cannot be compared. 

Similar tests dealing with the variabil ity on d ifferent time 
scales have been discussed by iMineshige et alJ lll994h and 
iMaccarone & Cor)rriM2002hl). 



4. Discussion 

4.1. The origin of the rms-flux relation 

lUttlev & McHardvl feOOlh have interprete d the rms-flux rela - 
tion in terms of the disk fluctuation model o flLvubarskiil(ll997l) . 
where variations in the accretion rate occur at various radii and 
propagate inwards to modulate the X-ray emission. The fact 
that variations at larger radii should have longer characteristic 
time scales naturally explains the fact that the short term rms 
variations are modulated by longer time scale flux changes. As 
the slower variations in the accretion flow propagate inwards 
they are able to mo dulate the variations on shorter time scales. 
iKotov etalJfcOOlh have extended this model by introducing an 
extended X-ray emitting region with a temperature gradient, to 
explain the frequency dependen t time lags between hard and 
soft bands observed in Cyg X-l JPottschmidt et all EoOO). 

Chur azov et alJ J200ll) have pointed out that in a standard 
thin disk, viscous damping will smooth accretion rate fluctua- 
tions on a viscous time scal e, i.e., before they reach the inner 
disk, implying that for the Lyubarskii i 19971) model to work, 
the fluctuations in accretion rate must happen in an optically 
thin or geometrically thick flow, such as an advection dom- 
inated accretion flow (ADAF) or another type of thick disk. 
We have shown that the linear rms-flux relation is seen in both 
high and low flux states, so that it very well might be associ- 
ated with some component of the system that is common to 
both states. Since the X-ray variability during the soft state 
appears to be associated with the power law emission rather 
than ther mal emission from the disk (which may be constant, 
IChurazov et allEpQlh . it seems likely that the rms-flux relation 
is associated directly with a corona which is present in both 
hard and soft states. We therefore speculate that the rms-flux 
relation in Cyg X-l (and probably also other objects showing 
linear rms-flux relations) may present evidence for accretion 
rate v ariations in a fluct uating coronal accretion flow (see, e.g., 
ISmith et all l2002t Ivan der Klisl 1200 lL and references therein). 
The suggestion that the corona itself accretes has been used to 
explain a number of ob servational fea tures of both AGN and 
galactic BH systems and since an optically 

thin coronal accretion flow is similar to an ADAF, such a model 
would not encounter the proble ms faced by a thin disk origin 
for the accretion variations ( Chura zov et 



by a thin 
all 1200 ll) . 



4.2. The long term rms-flux relation 

In the previous sections we have shown that the linear rms-flux 
relation seems to apply on very long time scales, since the av- 
erage 1-32 Hz rms of an observation scales linearly with the 
average flux of the observation. The hard and soft state data 
form two separate long term rms-flux relations, with different 
slopes &i ong (somewhat lower in the soft state, fc] ong =0.19 ver- 
sus fci ong = 0.22 in the hard state). Intriguingly, the flux-offsets 
Cion e m the long term rms-flux relations of both hard and soft 
states are very close to zero, and certainly of much smaller am- 
plitude than the typical offsets observed in the short term rms- 
flux relations of either state. The scatter in the hard state long 
term rms-flux relation is particularly small. 
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Combined with the almost-zero offset of the long term re- 
lation, this tight correlation implies that the fractional rms of 
the 1-32 Hz variability in the hard state is remarkably constant 
(similar to ki ong , i.e., ~22%). This result is rather surprising 
because the hard state PSD is clearly not stationary, i.e., the 
Lorentzians which successfully describe the PSD change sig- 
nificantly in both peak frequency and fractional rms (paper I). 
The most obvious example of such a change in the PSD can be 
seen in 1998 May (see paper I, Fig. 3) when the peak frequen- 
cies of the Lorentzians simultaneously increase, and the frac- 
tional rms of the L3 Lorentzian (which lies between 3-10 Hz) 
and the fractional rms of the total PSD both decrease signifi- 
cantly. However, despite such a marked change in the PSD, the 
data before and after 1998 May lie along the same long term 
hard state rms-flux relation. The fact that the fractional rms in 
the 1-32 Hz band seems undisturbed by significant changes in 
the PSD suggests some sort of conservation of fractional rms 
at high frequencies, even though the low-frequency (and hence 
the total) fractional rms changes significantly. Although the 
specific 1-32 Hz band chosen to measure the rms-flux relation 
is somewhat arbitrary, we point out that the bulk of the power 
in this band is associated with the L2 and L3 Lorentzians, and 
the band fully encompasses the range of frequencies between 
3-10 Hz in which a strong and well-defined anti-correlation be- 
tween Lorentzian peak frequency and fractional rms is clearly 
seen (see paper I, Fig. 6). We suggest, therefore, that some pro- 
cess which conserves fractional rms above a few Hz may be 
at work, with fractional rms in the L3 Lorentzian decreasing as 
the L2 Lorentzian moves further into the 1-32 Hz band, result- 
ing in a constant fractional rms in that band. We now examine 
how such a phenomenological model relates to the observed 
fundamental plane of short term rms-flux relation parameters 
k, C, and flux. 

4.3. The hard state k-C-(F) fundamental plane 

We have shown that the "fan shape" of the correlation between 
the rms-flux relation slope and its flux-offset in the hard state 
can be directly related to the long term rms-flux relation. The 
fact that k and C vary but the average 1-32 Hz rms must track 
the average flux explains the degeneracy in the correlation, 
which is the result of the various fluxes observed for a given 
pair of k and C values. Thus k, C, and (F) track out a funda- 
mental plane which describes the forms the rms-flux relation 
can take in the hard state. 

As noted by lUtrlev&McHardvl fcOOlh . C can be in- 
terpreted either simply as a constant-flux component to the 
lightcurve, or as a constant rms component. We concentrate on 
these two interpretations since they are the simplest possibili- 
ties, but it should be emphasized that there are numerous other 
forms the variability could take, e.g., a component of constant 
fractional rms. 

We have shown that, in the soft state at least, C can be neg- 
ative, which is only possible if the offset in this state represents 
a constant rms component. We argue that the existence of the k- 
C-{F) fundamental plane also suggests that C represents a con- 
stant rms component in the hard state. For example, consider an 



increase in the slope k between two observations, while keeping 
the average flux of the observations constant. The increase in 
k corresponds to an increase in the fractional rms of the linear 
rms-flux component of the lightcurve. In order to maintain the 
constant total fractional rms as given by the observed long term 
rms-flux relation, C must also increase. Therefore we suggest 
that both a constant rms component and a linear rms-flux com- 
ponent exist in the lightcurve and that the relative contribution 
of each varies in such a way that the sum of both components, 
the total fractional rms, is constant. 

As a constraint to possible models it should be mentioned, 
as already noted in Sect. 13.31 that the spectral shapes of k and 
C are generally similar, with C tending to a flatter/harder spec- 
trum, and there is no obvious correlation between spectral in- 
dex F and C, Ci ong , and k, respectively. 

4.4. Summary 

To summarize, in this paper we analyzed Cyg X-l lightcurves 
from 1996 to 2003 in an attempt to describe the long term evo- 
lution of X-ray rms variability. Our main results are 

1. The linear rms-flux relation is valid throughout all states 
and energy ranges observed. 

2. The slope of the linear rms-flux relation is steeper in the 
hard state than in the soft and intermediate states. 

3. The rms-flux relation is valid not only on short time scales 
(seconds) but also on longer time scales (weeks-months). 

4. In the hard state, the slope of the flux-rms-relationship is 
correlated with its flux intercept, but the correlation is in- 
creasingly degenerate at higher values of k and C, forming a 
fan-like distribution in the k-C plane. The hidden parameter 
governing the degeneracy is the long term average flux, (F), 
and the k-C-(F) relation maps out a "fundamental plane" in 
the low state. 

5. The intercept of the linear rms-flux relation is relatively 
constant and positive in the hard state, but is diverse in the 
soft state: in the 1996 soft state C is positive and slightly 
higher than during the hard state, in the atypical 2001/2002 
"soft state" C is decreasing relative to the hard state level, 
with several observations showing a negative C. During 
"failed state transitions", C is smaller than in the hard state, 
similar to the soft state behavior. 

6. The F-cr-T correlation can be confirmed for the hard state. 

7. The clear dichotomy of k in a hard and a soft state level in 
energy band 2 is weaker at higher energies. 
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Fig. A.l. In energy band 1 the rms-flux relation for observation 
2000 Nov 03, just as several other high flux observations during 
intermediate and soft states, displays an "arch-like" form. The 
given error bars are at the lcr level. 



Appendix A: Buffer overflows in RXTE binned 
data modes 

During phases of high soft X-ray flux, the rms-flux relation for 
energy band 1 (~2^4.6keV) shows a clear deviation from its 
typical linear behavior. These phases are characterized by a 
smaller rms than what would be expected when extrapolating 
the rms-flux relation measured for smaller fluxes (Fig. IA. li . 
We initially associated this breakdown of the relationship with 
changes i n the Comptonizing cor ona during phases of high lu- 
minosity (Gleissne r et all 12003). however, this interpretation 
was wrong. 

Instead, the "arches" are cause d by buffer over flows 
in th e PCA hardware (see also iGierlinski & Zdziarskii 
2003] Ivan Straaten et all l2003h . For the binned data mode 
B_2ms_8B_0_35_Q, with which most of our low energy data 
are taken, 4 bit wide counters are used to form the binned spec- 
trum during each 2 msec binning interval. During phases of 
very high flux, which are especially likely during short out- 
bursts possible during the failed state transitions, these buffers 
will overflow. Due to the softness of the X-ray spectrum dur- 
ing the failed state transitions, the likelihood for such over- 
flows is highest in the lowest energy band. The buffer over- 
flows become apparent when comparing the lightcurve of the 
binned data mode with the lightcurve as determined, e.g., from 
the GoodXenon data, which does not suffer from these prob- 
lems. Except for phases of high flux, these lightcurves agree. 
During phases of high flux, however, multiples of 16 events 
are missing in the lightcurve from the binned data. Since the 
buffer overflows only occur for a few bins of the high resolution 
lightcurve, the flux determination is barely affected, however, 
the rms is reduced and "arches" such as that shown in Fig. IA.ll 
are observed. For the computation of the rms-flux relationship, 
we therefore resorted to using energy band 2, where the overall 
photon flux is lower and buffer overflows do not occur. 



Note that GoodXenon data are not available for all obser- 
vations considered here, furthermore, the energy resolution of 
the GoodXenon data is not good enough to allow the study of 
other timing parameters. We note that the rms-flux relationship 
is especially sensitive to buffer overflows as the data are sorted 
according to flux. For the determination, e.g., of power spectra, 
longer lightcurves are chosen, and the buffer overflows are only 
apparent in a slight change of the rms level of the power spec- 
trum. Therefore, buffer overflows at the level seen in the data 
here only very slightly affect the normalization of the PSDs 
such as those used in paper I and do not have any influence on 
our earlier results. 
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